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Abstract 

We show that certain classical string configurations in AdS^ x T^'^ are chaotic. This 
answers the question of integrability of string on such backgrounds in the negative. We 
consider a string localized in the center of AdS^ that winds around two circles of T^'^. 
The corresponding dynamical system is equivalent to two coupled gravitational pendula 
and allows a very intuitive understanding. We find conclusive evidence of chaotic behavior 
by systematically analyzing the workings of the KAM theorem. We also show that the 
largest Lyapunov exponent is positive. 



1 Introduction 



One of the most interesting paradigms arising from the AdS / CFT [H [2l [3l H] correspondence 
has been the role of integrabihty (for an extensive review see |i5]). Integrabihty has modified 
our understanding of field theory to the point of reformulating, in some cases, the solution of a 
conformal field theory in terms of S-matrix elements of integrable systems. In practical terms, 
integrabihty has allowed exploration of sectors of the theory that would have been otherwise 
inaccessible with standard methods. 

In the context of the duality between strings on AdS^ x with Ramond-Ramond fluxes 
and = 4 supersymmetric Yang-Mills (SYM) with gauge group SU{N), integrabihty has two 
faces or two ways of appearing; they correspond to the two extreme values of the 't Hooft 
coupling. On the string theory side which admits a classical description at large 't Hooft 
coupling A, integrabihty of the classical sigma model on AdS^ x was suggested in p] and 
fully established with the inclusion of fermions in |7j. The field theoretic face of integrabihty 
originated out of the BMN limit, through the spin chain and relies on the form of the dilatation 
operator in an expansion around A — )■ 0. The key components of the integrabihty story are by 
now nicely summarized in [5]. 

More generally, the gauge/gravity correspondence provides a new method of attacking field 
theories at strong coupling. A natural hope is that the methods used in the standard duality 
could be developed to understand more realistic field theories and eventually apply them to 
theories in the same universality class of QCD. In this direction a first step was taken by 
Klebanov and Witten who considered a model with reduced supersymmetry. By taking the 
Maldacena limit of a stack of D3 branes at the tip of a conifold they arrived at a supergravity 
background of the form AdS^ x T^'^ [8]. This background is dual to an A/" = 1 superconformal 
field theory. 

The BMN limit of AdS^ x T^'^ was a clear indication that some of the properties of A/" = 4 
supersymmetric Yang-Mills could be found in certain sectors of superconformal A/" = 1 field 
theories with AdS dual [9l [ini HI] . More generally, it was shown in [10] that the corresponding 
BMN sector can be found in any gravity dual of the form AdS^ x where is a Sasaki- 
Einstein 5— manifold. Many classical string configurations found in AdS^ x have also been 
extended to analogous configurations in AdSr, x T^'^. For example, one can construct classical 
configurations of spinning strings in T^'^ [121 [13] • However, many unsuccessful attempts have 
been made to demonstrate that the classical string in AdS^ x T^'^ is integrable. In particular 
it can be established that the construction of [7] does not extend to AdS^ x T^'^ due to the 
fact that T^'^ is not a maximally symmetric coset space. The common consensus has been that 
some integrable structure is present but full integrabihty is, at least, elusive. 

Let us also remark that attacking the problem from the field theory side is rather hard as the 
corresponding field theory is defined at a strongly coupled conformal fixed point. This situation 
prevents a perturbative analysis from taking place. In other words, the fact that the fields have 
large anomalous dimensions in the conformal fixed point prevents a perturbative computation. 
An interesting study performed in [H] showed that the requirement of integrabihty of the full 
one loop dilatation operator in the scalar sector places very strong constraints on the field 
theory, so that the only soluble models correspond essentially to orbifolds of A/" = 4. 
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In this paper we show that in general the classical string in AdS^ x T^'^ can not be integrable. 
We settle this question not by our failure to show that the theory is integrable as has been 
the standard path so far. Rather, we present a simple classical configuration displaying chaotic 
behavior. The configuration in question is a string that wounds along two of the directions of 
T^'^; this configuration is inspired by the ring string in the Schwarzschild black hole in AdS^ 
considered in [15] but it is a lot simpler. As a result of the choice of Ansatz, the motion 
of the string reduces to the motion of two coupled oscillators that we analyze in detail. In 
a forthcoming work a similar albeit slightly more complicated construction is used to show 
non-integrability of string dynamics in AdS soliton background [T6]. 

The structure of the paper is as follows. In section|2]we present the Ansatz, its corresponding 
equations of motion and proceed to analyze various swaths of the phase space showing that 
the system presents some islands of integrability but is generically chaotic. To prove that the 
motion is chaotic we focus on the form of some Poincare sections; we use them to exemplify the 
workings of the Kolmogorov-Arnold-Moser (KAM) theorem. Namely, we show explicitly how 
the KAM tori get increasingly scattered as we change the energy. We also present a study of 
the the largest Lyapunov exponent in section |2j Section [3] contains a few comments about the 
field theory dual and section |4] contains our conclusions. 



2 Wrapped classical strings in T^'^ 

We start by considering the AdS^ x T^'^ background with metric 

ds"^ = (- cosh^ p dt^ + dp^ + sinh^ p dVtl 

1 ^ 1 ^ 

+ - ^{ddj + sin^ did4>1) + Q + XI 



i=l 1=1 

The Polyakov action for the string in the conformal gauge 
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S 



Ana' 



d'aGijdaX'd''X\ 



must be supplemented by the constraints 



drX' drX^ + d^X' d^X^ 



0, GijdrX' d,X^ = 0. 



(2.1^ 



(2.2) 



(2.3) 



We will consider solutions of the string localized in AdSr, at p = 0. The string will wrap the 
directions (pi and 02 in T^'^, a natural Ansatz will be 



(pi = aia, 02 = oi2o; 

t = t(r), tp = ^{T), e, = e,{r). 

Some of the coordinates involved have simple solutions t = E and ip 
equations are: 



(2.4) 

J, the nontrivial 
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— COS 6*2 I =0, 

ai \ 

— cos 6*1 1 =0, 

(2.5) 

The constraints yield 

= + ^ 5Z + ^ 5Z + ^ 5Z + cos e, COS (2.6) 

i=l i=l i=l 

It is easy to identify a natural classical analogy for our system. The above expression 
matches two gravitational pendula coupled through an interaction of the form cos^i cos02. We 
will exploit this mechanical analogy throughout our analysis. 

One of the most important results of XIX century mechanics was the identification of the 
role of the full phase space for characterizing a system, rather than the role of single trajectories. 
This is precisely the route that we will follow in the coming sections. The case for studying the 
full phase space has been made previously in the context of classical solution in the AdS/CFT 
correspondence [T7j. 



9i + «! sin 6*1 ( — cos 6i 
O2 + Ci2 sin 62 [ — cos 62 



2.1 Dynamics of the system 

In this section we discuss general aspects of the dynamics of the system. 

If one of the vanishes then the system reduces to a one dimensional model and hence is 
integrable in terms of Jacobi elliptic functions. Namely, for 02 = we have 62 = ar + h and 

+ ^ sin 201 = 0i(r) = am((r + C2)Jal/Q + C,, ). (2.7) 

The motion also becomes one dimensional if we start from sin{6i) = 0,6i = 0. Same holds for 

02. 

To investigate a more generic situation we look at the relevant potential term 

Vi9i, 62) = 4«? sin^ (0i) + a2 cos (0i) cos (^2) + ^aj sin" {62) . (2.8) 
io y io 

It can be easily seen from the eoms that only the ration the ratio — matters and both may 
be chosen to be positive. Local extrema for the above potential exist for (cti, a2 > 0), 



sin(0i) = O,sin(02) = (2.9) 
cos(0i) = 0, cos(02) = (2.10) 

sin(0i) = O,cos(02) = 2—, when |2— < 1| (2.11) 

cos(0i) = 2 — ,sin(0i) = 0, when |2— < 1| (2.12) 
a2 a2 
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Figure 1: Plot of the potential V{6i, 62) with ai = 1 and 02 = 2. 



sin ( 9i) 




Figure 2: Plot of sin(6'i(t)) with ai = 1 and 02 = 2 and a resting initial condition with 6*1(0) = 
0.2,6^2(0) = 0.1. The time evolution shows chaotic motion. 



Among these extrema the point 6*1 = 0, ^2 = is the local minimum (see Fig. Q. This 
suggests the existence of an stable fixed point near 9i = 0,02 = n. Near this point fluctuations 
of 61, 62 couple in quartic order and for sufficiently small value of those fluctuations the system 
behaves like a pair of decoupled harmonic oscillators. For such a motion we observe a nice 
periodic behaviour. For higher values of the energy the coupling between oscillators make the 
motion quasi-periodic and eventually leads to chaotic motion (Fig|2|. In the next section we 
study the details of the time evolution through phase space sections and discuss the issue of 
non-integrability. It is to be noted that for very high value of one angular momentum (i.e 6) the 
space virtually shrinks in those direction and the motion again becomes apparently integrable. 
This is related to the fact that for very high energy motion one may neglect the details of the 
potential. 

One interesting special case is ai = a2- Here we have a line of extrema along 61 = ±62, 
although this line of extrema always has an unstable direction. Any motion which starts 
with 9i = ±02 and 61 = ±62 remains on the line and essentially becomes a motion with one 
parameter and is integrable. Also the potential has enhanced symmetry. Near the minimum 
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of the potential there is an approximate rotational symmetry valid upto quartic order, which 
may lead to some interesting solutions. 



2.2 Poincare sections and KAM theorem 

In an integrable system the number of conserved quantities equals the number of degrees of 
freedom. A convenient way to understand those conserved charges is to look at the phase space. 
Let us assume that we have an integrable systems with N position variables qi with conjugate 
momenta pi. The phase space is 2N dimensional. Integrability means that there are conserved 
charges Qi = fi{p, q) which are constant of motion. One of them is energy. These charges define 
an A^-dimensional surface in the phase space which is a topological torus (KAM tori). The 2N 
dimensional phase space is nicely foliated by these N dimensional tori. In terms of action angle 
variables {Ii,6i) these tori just become surfaces of constant action. With each torus there are 
associated frequencies uji{li) which are the frequencies of motion in one particular action 
direction. 

What happens to those tori when the integrable hamiltonian is perturbed by a small non- 
integrable piece ? The KAM theorem states that most tori would be deformed but survive 
a small non-integrable deformation. However resonant tori which have rational ratio of fre- 
quencies, i.e. rrtiUi = with m E Q, will be destroyed. As the strength of the non-integrable 
interaction increases more tori gradually get destroyed. A nicely foliated picture of the phase 
space is no longer applicable and the time evolution may freely explore the whole phase space 
only constrained by energy. In such cases the motion becomes completely chaotic. 

To numerically investigate this gradual disappearance of foliation we look at the Poincare 
sections of our system. Our phase space has four variables 0i,92,pei,Pe2- If we fix energy we 
would be in a three dimensional subspace. Now if we start with some initial condition and 
time evolve, the motion would be confined to a two dimensional torus for the integrable case. 
This 2d torus intersects the sin(^2) = hyperplane at two circles. Taking gradual snapshots 
of the system as it crosses sin(6'2) = and plotting the value of (sin(6'i), cos(^i)p6)J, we can 
reconstruct those circles. Furthermore varying the initial boundary condition, in particular 
keeping (6'2(0) = vr,pe^ = 0) and varying 9i{0) to determine P6»2(0) from the energy constraint, 
we may hope to get a foliation structure typical for a integrable system. 

Indeed we see that for smaller values the of energy, a distinct foliation structure exists in the 
phase space [Fig 3(a)| . As we increase the energy some tori get gradually dissolved [Figs, 3(b) 



3(e)| . The tori which are destroyed sometimes broken down into smaller tori [Fig. 3(d)] eventu- 



ally the tori disappear and become a collection of scattered points known as cantori. However 
the breadths of these cantori are restricted by the undissolved tori and other dynamical ele- 
ments. Usually they do not span the whole phase space [Figs, 3(b) 3(d)| . For sufficiently large 



values of energy there are no well defined tori [Fig, 3(e)| . In this case phase space trajectories 



are all jumbled up and trajectories with very different initial conditions come arbitrary close 
to each other [Fig, 3(e)|. In terestingly for even higher energy order seems to have again form in 



the system [Figs ]3(f) 3(g)| . This is related to the fact that at very high energy we can neglect 



the potential. The mechanism is very similar to what happens in well known non-integrable 
periodic systems like double-pendulum models 
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(a) KAM tori 




1.0 -0.5 OO 



(b) Visible chaos 




-0.5 0.0 



(c) More into chaos 



-0.5 0.0 



(d) Even more into chaos 




(e) Completely chaotic 



" i. 



(f ) Formation of order 




(g) Ordered agian 
Figure 3: Poincare Sections. 
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2.3 Lyapunov exponent 



One of the trademark signatures of chaos is the sensitive dependence on initial conditions, which 
means that for any point X in the phase space, there is (at least) one point arbitrarily close 
to X that diverges from X. The separation between the two is also a function of the initial 
location and has the form AX{Xq, t). The Lyapunov exponent is a quantity that characterizes 
the rate of separation of such infinitesimally close trajectories. Formally it is defined as, 

A= hm hm -In^^^^ (2.13) 

r->ooAXo^or AX(Xo,0) ^ ^ 

In practice we use an algorithm by Sprott [201 EI] , which calculates A over short intervals and 
then takes a time average. We should expect to observe that, as time r is increased, A settles 
down to oscillate around a given value. For trajectories belonging to the KAM tori, A is zero, 
whereas it is expected to be non-zero for a chaotic orbit. We have verified such expectations 
for our case. We calculate A with various initial conditions and parameters. For apparently 
chaotic orbits we observe a nicely convergent positive A (Fig. |4]). 



3 Comments on the dual field theory 

The corresponding A/" = 1 superconformal gauge theory dual to AdS^ x T^'^ is known as the 
Klebanov-Witten theory and was originally described in [S|. The theory has flavor symmetry 
SU{2) X SU{2). The elementary degrees of freedom are denoted by the fields A and B, each a 
doublet of the factor groups SU{2) x SU{2) and with conformal anomalous dimension A^^^ = 
3/4. The gauge group is SU (N) x SU (N) and the two chiral multiplets A and B are respectively 
in the {N,N) and {N,N). The superpotential is 

W=^e'^e'''Tr[ABkAjBi], (3.1) 

with i = 1,2. 

The chiral operators analogue of the {X, Y, Z) operators in A/" = 4 SYM are in this case given 
by Tr{AB)^ with R-charge k and in the (|, |) representation of the flavor group SU{2) x SU (2). 

A discussion of the field content in relation to classical variables can be found, for example, 
in [9] whose notation we follow. As our starting point we can use a BMN-like operator of the 
form 



TrUM'] (3.2) 

This is an operator describing a spinning string in T^'^. In our case J corresponds to the angular 
momentum along the direction. This is a supersymmetric operator as it follows from the 
bounds of the supersymmetry algebra presented in [22] ■ Roughly speaking we associate angular 
momentum along 0i and 02 in T^'^ with the insertion of opertors of the form A and B. Let us 
schematically represent these operators as \ki, k,2 >■ These are, of course, not yet the kind of 
operators we are interested as we look for operators describing winding rather than momentum. 
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Figure 4: Motion of the string and the corresponding Lyapunov indices for a chaotic motion. We 
have chosen the initial condition 6*1(0) = 0.1, 6*2 = 0.1, 6*1 = 0,02 = 0.8. We see a convergent 
Lyapunov index with A ~ 0.245. 
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Thus we Fourier transform this operators along the same hues as for the magnon operatorn. 
Namely we consider operators of the form |ai,a2 >= X] l^i) ^2 >• To summarize, 

we take BMN-like operators and insert impurities that carry quantum numbers associated with 
angular momentum in the 0i and 02 directions and then perform a Fourier transform. 



4 Conclusions 

In this paper we have shown that certain classical configurations corresponding to a string 
winding along two of the angles of T^'^ exhibits chaotic motion for some range of initial condi- 
tions. This is enough to claim that the classical string in AdS^ x T^'^ can not be integrable and 
therefore we settle the important question of the existence of integrability beyond the known 
example of string theory on AdS^ x S^. We have thus showed that integrability seems to be 
a very peculiar property of S^, changing that manifold by what is considered the simplest 
next example (modulo orbifolds of S^) leads us into the realm of chaotic solutions. Given the 
simplicity of our dynamical system (two coupled gravitational pendula) we are very confident 
in the robustness of the numerical result which precisely matches what is expected from other 
text-book type of examples. Our example is actually very similar to the Henon-Heiles system 
and we exploited this similarity to make our case for chaotic behavior completely convincing 
following the standard literature on dynamical systems [TU [191 120] • 

As we saw in section |2] there are some islands, although of measure zero, of the phase space 
of the motion of the string in AdSr, x T^'^ where the system is integrable. For example, allowing 
the string to wrap only one direction in T^'^ yields an integrable system; it is only when we 
consider a configuration wrapping both directions (pi and 02 that we find chaotic behavior. 

It is worth remarking that the classical particle in ArfS's x T^'^ is integrable. The simplest 
way to see this is directly from the separability of the Laplacian. Of course, this is precisely 
the regime necessary for the study of supergravity modes. The implications of this situation 
- integrability of particles and chaos for strings - are translated in the dual field theory as 
corresponding to properties of operators dual to stringy states. Generically these are operators 
with very large quantum numbers and are not protected by symmetries. Recall, for example, 
that in the case of the BMN operators the classical string configuration consists of string shrunk 
to a point orbiting along a large circle in T^'^. The ground state operator is protected and the 
main point of BMN is that a subset of perturbations or insertions of operators remains semi- 
protected. What we anticipated due to the onset of chaotic behavior is that there are many 
operators for which insertions can not form a closed set. 

One interesting implication of our calculation is that it is easily generalizable to other Sasaki- 
Einstein manifolds such as Y^''^ and L^'^'''. As we have seen in section |2] a key property of the 
manifold that translated into the coupling of the corresponding oscillators was the existence of 
a nontrivial fibration. More precisely, the interaction contribution to the potential term which 



takes the form cos 6^1 cos 6*2 in equation (2.8) arises precisely because the space is a nontrivial 
f/(l) fibration over S"^ x S"^. The corresponding effect is present in the generic Sasaki-Einstein 
5-manifold which can be written as a nontrivial U{1) fibration over a Kahler-Einstein base. 
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